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ABSTRACT 

The detection of a nonthermal excess in the Coma cluster spectrum by two 
BeppoSAX observations analyzed with the XAS package (Fusco-Femiano et al. ) 
has been disavowed by an analysis (Rossetti & Molendi) performed with a differ- 
ent software package (SAXDAS) for the extraction of the spectrum. To resolve 
this discrepancy we reanalyze the PDS data considering the same software used 
by Rossetti & Molendi. A correct selection of the data and the exclusion of 
contaminating sources in the background determination show that also the SAX- 
DAS analysis reports a nonthermal excess with respect to the thermal emission at 
about the same confidence level of that obtained with the XAS package (~ 4. Sex). 
Besides, we report the lack of the systematic errors investigated by Rossetti & 
Molendi and Nevalainen et al. taking into account the whole sample of the PDS 
observations off the Galactic plane, as already shown in our data analysis of Abell 
2256 (Fusco-Femiano, Landi & Orlandini). All this eliminates any ambiguity and 
confirms the presence of a hard tail in the spectrum of the Coma cluster. 

Subject headings: cosmic microwave background — galaxies: clusters: individual 
(Coma) — magnetic fields — radiation mechanisms: non-thermal — X-rays: 
galaxies 



1. Introduction 

Nonthermal hard X-ray (HXR) emission was predicted in the seventies in clusters of 
galaxies showing extended radio regions, radio halos or relics, since the same radio syn- 
chrotron electrons can interact with the CMB photons to give inverse Compton (IC) X-ray 
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radiation (Perola & Reinhardt 1972; Rephaeli 1979). Several attempts to detect hard tails 
in the spectrum of a few clusters of galaxies were performed with various experiments (Baz- 
zano et al. 1984,90; Rephaeli, Grubcr & Rothschild 1987; Rephaeli & Gruber 1988; Rephaeli, 
Ulmer & Gruber 1994) that reported only upper limits to the nonthermal flux. A signifi- 
cant breakthrough in the measurement of nonthermal HXR emission was obtained thanks to 
the improved sensitivity and wide spectral capabilities of the BeppoSAX and Rossi X-Ray 
Timing Explorer {RXTE ) satellities. As pointed out by Petrosian (2003), the discovery of 
nonthermal HXR radiation has led to a remarkable increase of the theoretical investigations 
regarding the possible acceleration mechanisms and origin of the relativistic electrons re- 
sponsible for the nonthermal emission, although the presence of nonthermal phenomena in 
the intracluster medium (ICM) of some clusters was estabhshed decades ago (Willson 1970). 

Nonthermal HXR radiation was detected in excess of the thermal emission in the Coma 
cluster by a first BeppoSAX observation (Fusco-Femiano et al. 1999) using the Phoswich 
Detection System (PDS) and confirmed by a second independent observation with a time 
interval of about 3 yr (Fusco-Femiano et al. 2004; thereafter FF04). The presence of a second 
component in the X-ray spectrum of the cluster has been reported also by two RXTE obser- 
vations (Rephaeli, Gruber & Blanco 1999; Rephaeli & Gruber 2002). A2256 is the second 
cluster where a nonthermal excess has been measured by two BeppoSAX ohseiYations (Fusco- 
Femiano et al. 2000; Fusco-Femiano et al. 2005) and by RXTE (Rephaeh & Gruber 2003). 
At a lower confidence level, with respect to Coma and A2256, nonthermal HXR radiation has 
been detected by BeppoSAX in A754 (Fusco-Femiano et al. 2003). An upper hmit to the non- 
thermal fiux has been reported in A3667 (Fusco-Femiano et al. 2001), A119 (Fusco-Femiano 
et al. 2003) and A2163 (Feretti et al. 2001). For the last cluster a i?XT£' observation shows 
instead the presence of a HXR excess (Rephaeli, Gruber, & Arieli 2006). i?XTi? reports also 
some evidence of nonthermal emission by the Bullet Cluster (Petrosian, Madejski, & Lull 
2006). 

The PDS spectra of all the BeppoSAX observations were extracted using the XAS ver- 
sion 2.1 package (Chiappctti & Dal Fiumc 1997) specifically created to handle the PDS 
peculiarities. However, a PDS data analysis performed with a different software package 
(SAXDAS) does not report evidence for the presence of a hard tail in spectrum of the Coma 
cluster (Rossetti & Molendi 2004; thereafter RM04). 

To solve this contradiction we reanalyzed the PDS data of the two BeppoSAX ob- 
servations using the SAXDAS 2.0.2 package and in this letter we present our results and 
conclusions. In Sect. 2 we show the results of this analysis that are compared with those 
reported in FF04 using the XAS package and we examine the possible systematic errors 
investigated by RM04 and Nevalainen et al. (2004, thereafter NE04). In Sec. 3 we discuss 



-3- 



the reasons that have led to controversial results using different software packages. Finally, 
Sect. 4 is devoted to the conclusions regarding the presence of a hard tail in the spectrum 
of the Coma cluster. 

Throughout this Letter we assume a ACDM cosmology with Ho — 70 km s~^ Mpc~^ h-jQ, 
Q^rn — 0.3 and = 0.7. An angular distance of 1' corresponds to 29.0 kpc {zcoma — 0.0232). 
Quoted confidence intervals are at 90% level, if not otherwise specified. 

2. PDS Data Reduction and Results 

The Coma cluster was observed for the first time by the PDS instrument (Prontera 
et al. 1997) in December 1997 for ~91 ks and re-observed in December 2000 for ~300 ks. 
The pointing coordinates of BeppoSAX are at J(2000): a : 12'* 58"* 52^; 5: +27° 58' 54". 
To resolve the discrepancy between the results shown by FF04 and RM04 regarding the 
presence of a nonthermal component in the HXR Coma spectrum, we have re-analyzcd the 
PDS data with the SAXDAS package used by RM04. However, we start by reporting the 
procedure used in the XAS analysis and the main results obtained. More details can be 
found in FF04. 



2.1. XAS analysis 

The PDS data used in the XAS analysis were selected by a first automatic procedure 
followed by a visual check in order to eliminate all the remaining spikes that can affect signif- 
icantly the analysis results. The effective exposure times in the two PDS observations were 
44.5 ks and 122.2 ks, respectively (thereafter OBSl and 0BS2). Since the source is rather 
faint in the PDS band (~5 mCrab in 15-100 keV) a careful check of the background subtrac- 
tion was performed, making use of the default rocking law of the two PDS coUimators that 
samples ON/-FOFF, ON/-OFF fields for each coUimator with a dwell time of 96 sec (Pron- 
tera et al. 1997). When one collimator is pointing ON source, the other collimator is pointing 
toward one of the two OFF positions. Initially, we used the standard procedure to obtain 
PDS spectra (Dal Fiume et al. 1997); this procedure consists of extracting one accumulated 
spectrum for each unit for each collimator position. We then checked the two indepen- 
dently accumulated background spectra in the two different -I-/-OFF sky directions, offset 
by 210' with respect to the on-axis pointing direction (-I-OFF pointing: a : 12'* 58"* 57.8*; 
5: +24° 28' 55".l -OFF pointing: a : 12'* 58"* 47.0*; S : +31° 28' 54".7). The comparison 
between the two accumulated backgrounds (difference between the +OFF and -OFF count 
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Fig. 1. — PDS combined spectrum of the Coma cluster obtained with the XAS package 
(FF04). The continuous hne represents the best fit with a thermal component at the average 
cluster gas temperature of 8.11 keV (David et al. 1993). The error bars are quoted at the 
1(7 level. The spectrum starts at 15 keV. 
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rate spectra) showed that for OBSl the difference was compatible with zero (0.044 ± 0.047 
counts s~^ for a background level of 21.66 ± 0.02 counts s~^ in 15-100 keV), while for the 
longer, more sensitive 0BS2, there was an excess of 0.064 ± 0.021 counts s~^ (background 
16.76 ± 0.01 counts s~^ As reported in FF04, a careful check of possible variable sources 
in the PDS offset fields led the attention to the BL Lac source lES 1255+244, present in the 
+OFF field, that was observed by BeppoSAX on May 1998 in the framework of a spectral 
survey of BL Lacs by Beckmann et al. (2002). Because of the very short exposure time (~ 
3 ksec), our analysis of the source has determined only a 2a upper limit of 0.26 counts s~^ in 
15-100 keV, corresponding to 1.6 mCrab, however compatible with the background excess 
measured in 0BS2. Moreover, just in the center of the +OFF field is also present the 
extremely weak ROSAT source RX J125847. 1+242741. The presence of these contaminat- 
ing sources justified the decision to exclude the +OFF field in the background evaluation 
and consider only the -OFF field as the "uncontaminated" background for both the Coma 
observations. 

The combined spectrum, obtained by summing the spectra of the two observations 
(OBSl & 0BS2), shows a nonthermal hard excess with respect to the thermal component 
at the confidence level (c.l.) of ~ 4.8(T in the 20-80 keV energy range considering only the 
-OFF background spectrum. The spectrum is reported in Fig.l. FF04 also reported the 
confidence level of the nonthermal excess considering the standard procedure, i.e. by using 
as background the average of the spectra extracted from both the two offset fields. The 
c.l. is lower (~ 3.9cr) due to the contamination in the +OFF background. The average 
gas temperature used in the analysis is that measured by Ginga of 8.11+0.07 keV (David 
et al. 1993) in a field of view comparable to that of the PDS (FWHM ~ 1.3°). This value 
of the average temperature is confirmed by a determination of RXTE that reports a best-fit 
temperature of 7.90+0.03 keV (Rephaeli & Gruber 2002) in a field of view of ~ 1° comparable 
to the field of view of Ginga and PDS. 



2.2. SAXDAS analysis 

The first combined spectrum was extracted considering the PDS data resulting by the 
automatic selection operated by the SAXDAS package %. The confidence level of the excess, 
in the range 20-80 keV and for an average gas temperature of 8.11 keV, is ~ 2.9cr (observed 



^The ~20% variation in the PDS background is due to the BeppoSAX orbital decay: the lower orbit for 
0BS2 increased the shielding to ambient particles, therefore lowering the diffuse background. 
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count rate = 0.1717±0.0146 counts s"^ , model predicted rate = 0.1295 counts s~^ ) 
taking into account only the uncontaminated -OFF position. The effective exposure time 
is 42.8 ks in OBSl and 119.3 ks in 0BS2 for a total exposure time of 162.1 ks. To resolve 
the discrepancy between these results and those obtained by FF04 with a XAS analysis 
(see previous section), a new combined spectrum was extracted considering the same time 
windows used in the XAS analysis. These time windows are different from those used by 
RM04 in their SAXDAS analysis. The reason is in the different criterion for eliminating the 
spurious spikes produced by the charged particles that interact with the PDS detectors: a 
semi-automated procedure (i.e. an automated procedure followed by a visual check) for XAS 
and a completely automated procedure for SAXDAS (Guainazzi & Matteuzzi 1997). With 
the same time windows used by XAS, the c.l. of the excess in the SAXDAS spectrum raises 
at ~ 4.2cr (observed count rate = 0.1902±0.0148 counts s^^ , model predicted rate — 0.1280 
counts s~^ ), while the total exposure time reduces to 160.9 ks (it was 162.1 ks). 

The two packages, starting from the same time windows, produce a difference in the 
total effective exposure time of 5.8 ks {txAs = 166.7 ks, tsAXDAs = 160.9) essentially due 
to: a) the different earth angle above which data are selected (the earth angle is defined as 
the angle between the earth limb and the BeppoSAX '^omiing direction). The earth angle 
is 5° for XAS and 10° for SAXDAS; b) the removal of observational time after any South 
Atlantic Geomagnetic Anomaly (SAGA) passage. The SAXDAS package removes 5 minutes 
after any passage, while XAS ehminates the time necessary to reach the correct voltage of 
the instruments after their shut-down during the SAGA passage. 

In conclusion, we have a lower exposure time for the SAXDAS analysis which may imply 
a lower c.l. of the nonthermal excess. To quantify this difference we have extracted a further 
combined spectrum using the same time windows selected in the XAS analysis but imposing 
an earth angle of 5° in the SAXDAS package. The total effective time exposure of the two 
observations is now tsAXDAS — 169.1 ks instead of 160.9 ks. The excess results to be at the 
c.l. of ~ 4.6cr (observed count rate = 0.1944±0.0144 counts s~^ , model predicted rate = 
0.1280 counts s~^ ). This HXR spectrum is reported in Fig. 2. The difference in the initial 
energy of the SAXDAS and XAS spectra (SAXDAS spectrum starts at 12 keV, while XAS 
at 15 keV) gives insignificant variations in the c.l. values of the nonthermal excess. 



2.3. Possible systematic errors 

The systematic errors examined by RM04 and NE04 have already been discussed by 
us in the PDS analysis of A2256 (Fusco-Femiano, Landi & Orlandini 2005), but considering 
their importance in the analysis results we intend to report here again the main parts of the 
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Fig. 2. — PDS combined spectrum of the Coma cluster obtained with the SAXDAS package 
(see text). The continuous hne represents the best fit with a thermal component at the same 
average cluster gas temperature of Fig. 1. The error bars are quoted at the la level. The 
spectrum starts at 12 keV. 
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discussion. 

1) RM04 claim the presence of an "Instrumental background residual" (see sect. 2.1 of 
their paper) derived from the analysis of 15 "blank fields", i.e. fields which do not contain 
sources showing significant emission in the PDS energy range. By summing the spectra 
from these observations they find that the spectrum differs from zero: the count rate is 
(1.45 ±0.77) X 10~^ counts in the 12-100 keV energy range. This seems to indicate that 
the background in the ON position is larger than that in the ±OFF positions producing 
an instrumental contribution not removed by the background subtraction procedure. This 
effect has been studied by Landi et al. (2005) considering the complete sample of 868 PDS 
pointings with galactic latitude \b\ > 15°, and selecting the 15-100 keV net count spectra for 
which there is source detection below la (that is, "blank fields"). These spectra have been 
summed imposing a net exposure greater than 20 ksec. A net count rate of (1.67±5.30) x 10~^ 
counts has been derived, consistent with the definition of "blank field" . Also NE04 do 
not report evidence for an instrumental residual. 

2) The other effect evaluated by RM04 regards the systematic differences between the 
background fields. They analyze a sample of 69 observations whose target is outside the 
galactic plane and with a long exposure time (see Appendix of their paper) . RM04 find that 

the mean value of the difference between ON and the two -OFF and +OFF sky positions is 
significantly different from zero and positive. Also this effect has been investigated by Landi 
et al. (2005) on the whole sample of PDS observations. The obtained value of (5.3±6.3) x 10~^ 
counts s~^ is consistent with no contamination at all. We presume that the value found by 
RM04 could be due to the small sample of observations they considered. NE04, analyzing 
a larger sample of data with respect to that used by RM04 (164 PDS observations), found 
a systematic difference between ON and the two offset pointings that cancels out in the 
standard usage of both offsets. 

3) In addition, NE04 introduce a systematic error in the net source count rates due 
to unresolved and not significantly detected point sources present in the PDS field of view. 
They find that an excess of 0.019 counts s~^ has to be added to the net count rate spectra (no 
errors are given on this measurement) , when the standard method of background evaluation 
is used, and 0.027 counts s~^ has to be added when the background is evaluated from only 
one offset field. We performed this same analysis on a set of 868 observations (NE04 use 164 
fields) and find that the contribution of background fiuctuations due to unresolved and not 
significantly detected sources in the offset fields (in other words, the PDS confusion limit) 
introduces a variance crji^^ = (9.5 ± 10.3) x 10"'' (counts s~^)^ consistent with zero (Fusco- 
Femiano, Landi & Orlandini 2005). Therefore PDS data are not affected by this systematic 
effect. 
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3. Discussion 

Fig. 2 shows that also the PDS combined spectrum extracted with the SAXDAS pack- 
age evidences the presence of a nonthermal excess with respect to the thermal emission at 
about the same c.l. reported by the XAS analysis (FF04) (~ 4.8a" for XAS and ~ 4.6a" 
for SAXDAS). The discrepancy between the results reported in FF04 and RM04 is mainly 
due: a) to an accurate selection of the events. In fact, we have shown that the SAXDAS 
analysis of the same PDS time windows used with the XAS software leads to a significant 
increase of the c.l. of the excess (from ~ 2.9a" to ~ 4.2a); b) to a correct determination of 
the background. The XAS and SAXDAS spectra shown in Fig.l and Fig. 2, respectively, 
are obtained considering only the -OFF background. The exclusion of the -l-OFF position in 
the background determination implies a more pronounced detection of nonthermal radiation 
from the Coma cluster. The check regarding the presence of contaminating sources in the 
+OFF field was not operated by RM04. 

An increase of ~ 0.4a in the c.l. value is obtained adopting in the SAXDAS analysis 
the same earth angle of 5° used in the XAS procedure. The lower earth angle implies a 
greater exposure time that justify an improvement in the c.l. value of the nonthermal excess. 
Considering in the SAXDAS analysis the same time windows and the same earth angle used 
in the XAS analysis the difference in the total exposure time is 2.4 ks in favour of the 
SAXDAS package (txAs = 166.7 ks, tsAXDAs = 169.1 ks) due to the different removal of 
the time operated by the two packages after any SAGA passage. The XAS package results 
to be more conservative waiting more time than SAXDAS to allow the PDS high voltage to 
reach the correct levels. So, some spurious events could be present in the SAXDAS analysis, 
in particular for the longer 0BS2 observation. 

In order to reproduce the results obtained by RM04 we have computed the c.l. of the 
excess for OBSl and 0BS2 spectra (RM04 do not report the c.l. value for the combined 
spectrum) using the automatic selection of the events operated by the SAXDAS package and 
the standard procedure for the background determination. The c.l. values are 2.90a and 
1.34a" for OBSl and 0BS2, respectively, in the energy range 25-80 keV and for an average 
gas temperature of 8.21 keV. These values are not much distant from those reported in Table 
2 of RM04 (2.84a for OBSl and 1.11a for 0BS2)@ confirming that the selection procedure of 
the PDS events and the exclusion of the +OFF position for the background determination 
are the main reasons of the different results reported in FF04 and RM04. The PDS data 
analysis of very weak sources like Coma in the HXR band requires a rigorous selection of 



^These values of the c.l. of the excess refer to the line "No Subtraction" (no subtraction of the instrumental 
residual that Landi et al. (2005) find consistent with zero; see Sect. 2.3). 
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the events in order to eliminate the presence of any spikes able to introduce noise that hides 
the presence of a nonthermal excess with respect to the thermal radiation. The presence of 
contaminating sources in the offset fields does not allow the use of the standard procedure 
for the background evaluation. 

4. Conclusions 

We have shown that the presence of a nonthermal excess with respect to the thermal 
emission in the spectrum of the Coma cluster does not depend on the used software pack- 
age (XAS or SAXDAS) for the PDS data analysis. The spectra extracted with XAS and 
SAXDAS, reported in Fig. 1 and 2, respectively, show both a nonthermal excess at about 
the same c.l. value when in the SAXDAS analysis the same time windows used in the XAS 
analysis (FF04) are adopted and the +OFF sky direction is not taken into account in the 
background determination for the presence of contaminating sources. The systematic ef- 
fects claimed by RM04 and NE04 can be excluded considering the whole sample of the PDS 
observations off the Galactic plane. 

This re-analysis of the PDS data, using the SAXDAS package, explains the different 
results reported in FF04 and RM04 confirming the presence of a nonthermal component in 
the Coma cluster spectrum. 
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